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We find that the amplitude of the R xx radiation-induced magnetoresistance oscillations in 
GaAs/AlGaAs system grows nonlinearly as A oc P a where A is the amplitude and the exponent 
a < 1. This striking result can be explained with the radiation-driven electron orbits model, which 
suggests that the amplitude of resistance oscillations depends linearly on the radiation electric field, 
and therefore on the square root of the power, P. We also study how this sub-linear power law 
varies with lattice temperature and radiation frequency. 
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I. INTRODUCTION 

Superconductivity^] and quantum Hall effects[2, 3] are 
known to present two distinct and interesting examples 
of zero-resistance states in material systems. Another 
interesting example occurs in the quantum Hall two- 
dimensional electron system (2DES) when the 2DES is 
irradiated with microwave and terahertz band radiation 
in the presence of a transverse magnetic field. Here, 
it becomes possible to photo-excite into zero-resistance 
states in this low dimensional system. [4, 5] Photo-excited 
transport in the 2DES has become a topic of experi- 
mental and theoretical interest since this observation. [4- 
44] Periodic in B~ x radiation- induced magnetoresistance 
oscillations, which lead into the radiation-induced zero- 
resistance states, are now understood to be a consequence 
of radiation-frequency (/) and magnetic field (B) depen- 
dent, scattering at impurities [24-28] and/or a change 
in the distribution function[6, 31, 38]. And, vanishing 
resistance at the oscillatory minima is explained as a fea- 
ture of negative resistance instability and current domain 
formation. [25, 34, 42] In spite of the progress, there re- 
main many aspects that could be better understood from 
including, for example, the growth of the oscillatory effect 
vs. the radiation intensity, P. Here, a number of works 
have numerically evaluated the radiation-induced mag- 
netoresistance oscillations for several P and graphically 
presented the results, [24, 32, 35] while Dmitricv ct al.,[31] 
have suggested that the correction to the dark dc con- 
ductivity is linear in P. [31] A comparison of experiment 
with theory, so far as the P-dependence is concerned, 
could help to identify the importance of the invoked- 
mechanisms in these theories. [36] 

Thus, we examine the growth of the radiation- induced 
magneto-resistance oscillations with P. Experiment in- 
dicates that the amplitude A of the radiation-induced 
oscillatory diagonal resistance (R xx ), grows nonlinearly 



with P and can be described by A oc P a where the expo- 
nent a < 1. At the same time, according to experiment, 
a also depends on the lattice temperature T and radi- 
ation frequency /. Since such non-linear growth of A 
with P, T, and / had not been predicted, we propose 
a theoretical explanation. The explanation utilizes the 
radiation-driven electron orbit model,[32] where radiation 
forces the orbit center of the Landau states to move back 
and forth in the direction of the radiation electric field 
at the frequency of radiation, and the R xx oscillations 
reflect the periodic motion of the electron orbits center. 
This theory establishes that the amplitude {Ai, an ) of the 
orbit center motion, sets the amplitude A of the R xx os- 
cillations, and that AL an is proportional to the radiation 
electric field. Since the radiation power P is proportional 
to the square of the radiation electric field, it follows that 
R xx and A will depend on the square root of P, at the 
lowest temperatures. The experiments indicate that a is 
close to 0.5 at lower temperatures. 



II. THEORETICAL MODEL 

The above mentioned theory[32, 33, 39] was proposed 
to explain the R xx of an irradiated 2DES at low B. We 
obtained the exact solution of the corresponding elec- 
tronic wave function[32, 39, 45-48]: 



OC (f> n (x - X - Xd(t),t) 
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, where <j> n is the solution for the Schrodingcr equation 
of the unforced quantum harmonic oscillator, X is the 
center of the orbit for the electron motion, x c i(t) is the 
classical solution of a forced harmonic oscillator: 
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where e is the electron charge, 7 is a phenomenologically- 
introduccd damping factor for the electronic interaction 
with acoustic phonons, w c the cyclotron frequency, E 
the radiation electric field. Then, the obtained wave 
function is the same as the standard harmonic oscilla- 
tor where the center is displaced by x c i(t). Thus, the 
orbit centers are not fixed, but they oscillate harmoni- 
cally at the radiation field frequency w = 2nf. Then, 
by this micro(MW)-driven periodic motion, electrons in 
their orbits (quantum oscillators) interact with the lat- 
tice ions being damped and emitting acoustic phonons. 
In the x c i expression, 7 represents this damping. 

This radiation — driven behavior affects dramatically 
the charged impurity scattering and eventually the con- 
ductivity. Then, first we calculate the impurity scatter- 
ing rate Wn,m between two oscillating Landau states 
^n, and <3/m[32, 47, 49, 50]. Next we find the average 
effective distance advanced by the electron in every scat- 



tering jump: AX 
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A-Lan cos idt [32, 39, 47], where 



t = 1/Wn,m is the scattering time. Finally the longitu- 
dinal conductivity a xx is given by: 
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being E the energy. To obtain R xx we use the relation 
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and the amplitude of resistance oscillations depends lin- 
early on the radiation electric field. 

P can be related with R xx through the well-known 
formula that gives radiation intensity / (power divided 
by surface) in terms of the radiation electric field Eq: 
I = t}C€qEq, where c is the speed of light in vacuum and 
6q is the permittivity in vacuum. If we want to express 
only the power in terms of the radiation electric field 
we have to take into account the sample surface. In the 
particular case of GaAs if S is the sample surface, we can 
readily obtain: 



P — -^ c GaAs(-GaAs^E S 



(5) 



where cq cl a s is the speed of light in GaAs and tGaAs is 
the GaAs dielectric constant. Accordingly, Eq oc \f~P . 
Then, substituting in the expression of Alcm, we obtain 
that R xx varies with P following an square root law: 
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Thus, the R xx response will grow as the square root 
of P, i.e., A (x P 5 . From the expression of R xx we 
can study the influence of w and T on the sublinear rela- 
tion. In the case of w, the calculation is straightforward 



(see expressions of R xx or AL an ). In the case of T, wc 
introduce a microscopic model, which allow us to ob- 
tain the damping parameter 7 and its dependence on T 
[32, 33, 50]. According to the model, 7 is proportional 
to the scattering rate of electron-acoustic phonon inter- 
action, being eventually linear with T[50]. 



EXPERIMENTAL SET UP AND RESULTS 

Low frequency lock-in based electrical measurements 
were carried out at T < 1.5K with the samples immersed 
in pumped liquid-helium and mounted near the open end 
of a microwave waveguide. [4, 13] The W (Wegscheider)- 
GaAs/AlGaAs single heterostructures were nominally 
characterized by an electron density, n = 2.4 x 10 n cm -2 
and a mobility of \i = 10 7 cm 2 /Vs. Results are reported 
here for measurements on a set of three W-specimcns. 

To examine the growth of the radiation-induced os- 
cillations with P, Fig. I presents the AR XX of specimen 
Wl for several P at 50GHz. Also shown in the figure are 
data-fits to A.Rf% = -Aexp{-\/B)sin(2irF/B). Here, 
a slowly varying background, approximately equaling the 
dark trace, was removed from the photo-excited R xx data 
to realize Ai? a;x .[29] Although this fit function includes 
three parameters, A, A, and F, the oscillation period in 
B^ 1 is independent of the radiation-intensity, and thus, 
F is a constant. Further, the damping constant, A, [9] 
turns out to be insensitive to P, see Fig. 1(b). Thus, the 
main free parameter in the fit-function is the amplitude, 
A, of the oscillations. In Fig. 1(c), we exhibit the fit 
extracted A vs. P for Wl and W2. The figure shows a 
sub-linear growth of A with P. Also shown are power law 
fits, A = A P a . Here, a = 0.63 and a = 0.64 for Wl and 
W2, respectively, at 1.5K and 50GHz. Aq varies between 
Wl and W2 because the effective intensity attenuation 
factor varies between the two experiments. 

For a third W-specimen labeled W3, Fig. 1(d) and 
1(e) report the influence of the temperature on the 
growth of A vs. P, where A is extracted, as before, 
from line-shape fits of the oscillatory data to AR^ = 
-Aexp{-\/B)sin{2TiF/B). Here, Fig. 1(d) and 1(c) 
indicate that, as expected, at a constant P, A grows 
with decreasing T both at / = 50GHz [Fig. 1(d)] and 
/ = 46GHz [Fig 1(e)]. Further, the figures show that 
the A vs. P curves exhibit greater curvature at lower 
temperatures. The A vs. P have been fit once again to 
A = AoP a ; the fit-extracted a and Aq have been sum- 
marized in tabular form within Fig. 1(d) and Fig. 1(e). 
These fit-extracted a have also been plotted vs. T in 
the inset to these figures. These insets suggest that a 
decreases with decreasing temperatures, consistent with 
the observed increased non-linearity at lower tempera- 
tures. Note that, at / = 50GHz, all three W-specimens, 
exhibit the same a, within uncertainties, at the lowest 
pumped 4 He temperatures. In addition, a comparison 
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of the a reported in Fig. 1(d) and Fig. 1(c) also sug- 
gests that reducing the microwave frequency / at a fixed 
T tends to reduce the a, i.e., increase the non-linearity. 
Thus, the experimental results presented here suggest a 
nonlinear power law in a regime characterized by mod- 
est excitation. This peculiar behavior can be theoret- 
ically explained and the experimental results recovered 
by means of the radiation-driven electron orbits model. 
Next we present the calculated results that seem to rea- 
sonably agree with the present experiments. 

CALCULATED RESULTS 

Fig. 2 presents the calculated R xx as a function of B 
over < P < 3.7mW at / = 50GHz and T = IK. As 
in experiments, we obtain multiple oscillations for all P 
except, as expected, for P — 0. It is clear the progressive 
collapse of R xx oscillations as P decreases. The explana- 
tion is straightforward according to our theory and it is 
given by the expression of R xx . Thus, for a decreasing 
(increasing) P the R xx response is progressively smaller 
(larger). The key issue here, and the challenge for the dif- 
ferent available theories, is to deduce the rate of increase 
(decrease) of R xx oscillations for increasing (decreasing) 
P. This is shown in the inset of Fig. 2, where we present 
the calculated amplitude versus P. As in experiments, we 
obtain a clear sublincar growth of A with P: A oc P a . 
We would expect an exponent a around 0.5 since the 
model establishes an square root dependence. Thus, we 
have performed a fit of the calculated values, obtaining 
A = 1.66P 5 , very close to a square root law. 

The dependence of the sublincar law on T and fre- 
quency is presented in Fig. 3. In Fig. 3a we represent 
the calculated amplitude A of the main R xx peak versus 
P for a frequency of 50GHz and T values: T = 0.75 A', 
1.0A, 1.5K and 2. OX'. We obtain a qualitatively simi- 
lar behavior as in experiments. Thus, we observe that 
for a constant P, A increases (decreases) with decreasing 
(increasing) T. The physical explanation is a follows. If 
one increases T, the interaction between electrons and 
lattice ions also increases giving rise to a more intense 
damping. This implies progressively smaller amplitudes 
and R xx oscillations tend to vanish. On the other hand, 
it is very clear that for all T used in the calculation, the 
dependence between A and P is sub-linear and that, at 
least qualitatively, the curves show a bigger curvature for 
decreasing T. To check this curvature, we have made a 
fit to the calculated values for each T observing that the 
exponent of the fits slightly increases for increasing T. 
In all cases, the exponent is always close to 0.5. This 
outcome is expected since, strictly speaking, our model 
does not predict any variation of the exponent with T; 
our theoretical framework imposes a square root relation. 
Yet, experimental results show a much faster increase of 
the exponent with increasing T regarding the calculated 



ones. The latter show a much slower variation with T, 
although in the same direction. At present, this model 
can not explain this quantitative discrepancy of the vari- 
ation of a with T. This discrepancy could be the result 
of an statistical effect of the fit or the reflex of a real 
physical process. In the latter case, this would imply an 
extension of the current theory that should start from 
the square root law obtained here. We consider that the 
square root dependence is a solid result reflecting a rea- 
sonably good agreement between theory and experiments 
a low temperatures. 

A similar behavior is presented in Fig. 3b where we 
show the calculated values of A versus P for different 
frequencies. We observe that, as in experiments, A de- 
creases as the the frequency increases. The explanation 
comes from the expression of R xx where w shows up in 
the denominator. Then, for a constant P, a larger w 
gives rise to smaller amplitudes. We have fit these curves 
and, as expected, the exponent of the fit slightly increases 
for increasing frequency whereas the pre-factor decreases. 
Yet, experiments show a much faster increase. This 
discrepancy is similar to the case of the T-depcndcncc. 
Thus, the explanation given in Fig. 3a, applies identi- 
cally here too. 

All calculated results presented in Figs. 2 and 3, have 
been obtained directly through equation (6). 



CONCLUSIONS 

In summary, experimental results indicate a nonlinear 
growth in the amplitude A of radiation-induced magneto- 
resistance oscillations with P. These results can be ex- 
plained with the radiation-driven electron orbits model. 
According to this model, the amplitude of the radiation- 
induced oscillations should be proportional to the square 
root of radiation power, which implies a sub-linear rela- 
tion. We have also studied how this sub-linear law varies 
with lattice temperature and radiation frequency obtain- 
ing only a qualitative good agreement between theory 
and experiments. The obtained quantitative discrepan- 
cies can not be explained by our theory. The origin of 
them could be an statistical effect or a real physical ef- 
fect. In the latter case, this would imply an extension of 
the present theoretical model, always starting from the 
square root law here obtained. 
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FIG. 1: (a) For Wl, AR XX is exhibited at / = 50GHz. Also 
shown are fits to an exponentially damped sinusoid, (b) A is 
plotted vs. P for Wl and W2. (c) The lineshape amplitude, 
A, is plotted vs. P for Wl and W2. Also shown are fits, 
A = A P a , which suggest a = 0.63 and a = 0.64 for Wl 
and W2, respectively. (d) At / = 50GHz, the amplitude, A, 
is plotted vs. P for T = 1.4ff, 1.8 A* and 2.2A" for W3. Also 
shown are fits to A = yloP Q . The fit-extracted a and 
are presented are presented in tabular form within the figure. 
The inset shows the variation of a with T at 50GHz. (e) 
At / = 46GHz, the amplitude, A, is plotted vs. P for T = 
1.4A, 1.75A and 2.6A for W3. The inset shows the variation 
of a with T at AGGHz. 
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FIG. 2: Calculated R xx vs. B at / = 50GHz, T = Iff, 
and < P < 3.7mW. We obtain a progressive quenching 
of R xx oscillations with decreasing P, leading into the dark 
curve (P — OmW). The inset shows the calculated amplitude 
of the main R xx peak of Fig. 2 (B = 0.15T) versus P. A 
increases sub-linearly with P. The line corresponds to a fit of 
the calculated values with A = 1.66P ' 55 which, as expected, 
is very close to a square root law. 
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FIG. 3: (a) Calculated amplitude A of the main R xx peak 
versus P for / = 50GHz and the T values: T = 0.75K, 1.0K, 
1.5K and 2.0K. Note that the dependence between A and P 
is sub-linear for all T. Also shown are the fits performed for 
each T. Observe that for increasing T, the exponent of the fit 
increases as the pre-factor decreases, (b) Calculated A versus 
P for radiation frequencies, 40GHz, 50GHz and 70Ghz and 
T — IK. We observe that A decreases as / gets larger. A fit 
of the A vs. P curves indicates that the exponent increases 
with increasing frequency, whereas the pre-factor decreases. 



